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Supersymmetric extensions of the standard model of particle physics assuming the gravitino to 
be the lightest supersymmetric particle (LSP), and with the next-to-LSP decaying to the gravitino 
during Big Bang nucleosynthesis, are analyzed. Particular emphasis is laid on their potential to solve 
the " 7 Li problem" , an apparent factor 2 — 4 overproduction of 7 Li in standard Big Bang nucleosyn- 
thesis (BBN), their production of cosmologically important amounts of 6 Li, as well as the resulting 
gravitino dark matter densities in these models. The study includes several improvements compared 
to prior studies concerning NLSP hadronic branching ratios, NLSP dark matter densities, the eval- 
uation of hadronic NLSP decays on BBN, BBN catalytic effects, as well as updated nuclear reaction 
rates. Heavy gravitinos in the constrained minimal supersymmetric standard model (CMMSM) are 
reanalyzed, whereas light gravitinos in gauge-mediated supersymmetry breaking scenarios (GMSB) 
are studied for the first time. It is confirmed that decays of NLSP staus to heavy gravitinos, while 
producing all the dark matter, may at the same time resolve the 7 Li problem. For NLSP decay 
times ss 10 3 sec, such scenarios also lead to cosmologically important 6 Li (and possibly 9 Be) abun- 
dances. However, as such scenarios require heavy ^ 1 TeV staus they are likely not testable at the 
LHC. It is found that decays of NLSP staus to light gravitinos may lead to significant 6 Li (and 9 Be) 
abundances, whereas NLSP neutralinos decaying into light gravitinos may solve the 7 Li problem. 
Though both scenarios are testable at the LHC they may not lead to the production of the bulk of 
the dark matter. A section of the paper outlines particle properties required to significantly reduce 
the 7 Li abundance, and/or enhance the 6 Li (and possibly 9 Be) abundances, by the decay of an 
arbitrary relic particle. 



I. INTRODUCTION 

Supersymmetry is probably the best studied extension 
of the standard model (SM) of particle physics. Softly 
broken supersymmetry seems attractive as it not only 
may solve the hierarchy problem but also lead to gauge 
coupling unification at the GUT scale and electroweak 
symmetry breaking by radiative corrections. Addition- 
ally, supersymmetry predicts the existence of new parti- 
cles whose lightest (LSP), if stable due to some R-parity, 
may be naturally produced in the right abundance to 
provide the cosmological dark matter. Here the case of 
neutralino LSPs has been widely discussed over at least 
one-and-half decades, whereas the case of gravitino LSPs 
has received the same widespread detailed attention only 
since about five years 0, 0, i, S S S, 0, S I, M EL El 
Gravitino dark matter is special due to 
its superweak interactions with ordinary matter. It may 
be produced at fairly high temperatures 10 5 — 10 s GeV 
during reheating after inflation and, in case it is the LSP, 
by the decay of the next-to-lightest supersymmetric par- 
ticle (NLSP) to its SM partner and a gravitino. Due to 
the gravitationally supressed interactions of the gravitino 
the decay with decay time 
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typically takes place during or after Big Bang nucle- 
osynthesis (BBN). Here Mnlsp and m G are NLSP- and 



gravitino- mass, respectively, and M p \ denotes the re- 
duced Planck mass, [17]- Gravitino dark matter is there- 
fore subject to constraint, as it may leave its imprint in 
the light element abundance yields. 



Prediction of the light element abundance yields of 
2 H (D), 3 He, 4 He, and 7 Li by an epoch of Big Bang 
nucleosynthesis is one of the big successes of the hot 
Big Bang model. It has not only lead to the realiza- 
tion that the Universe has expanded by at least a fac- 
tor 10 10 in its past, but also that the bulk of the dark 
matter must be non-baryonic. In its standard version 
BBN is reduced to a model of one parameter only, the 
fractional contribution of the baryonic density to the 
critical one, Sl^h 2 (where h is the present Hubble con- 
stant in units of 100 km s _1 Mpc _1 ). Five years of ob- 
servations of the cosmic microwave background radia- 
tion (CMBR) anisotropies by the WMAP satellite [H| 
have lead to the comparatively accurate determination 
of n b h 2 fa 0.02273 ± 0.00062 (WMAP-only). This in- 
dependent estimate of Vl^ti 2 has promoted a comparison 
of observationally inferred primordial light clement abun- 
dances with those predicted in standard BBN to an inde- 
pendent cross-check of the cosmic SM. Such a comparison 
is very favorable for D, broadly consistent yet somewhat 
inconclusive for 3 He and 4 He, and significantly discrepant 
for the isotope of 7 Li. In particular, the primordial 7 Li/H 
ratio is commonly inferred from 7 Li absorption lines in 
the atmospheres of low-metallicity halo, or globular clus- 
ter stars. Most determinations yield 7 Li/H ratios in the 
range 1 - 2 x 10" 10 [H, M, HE HI HI , essentially the 
same value as that determined by the Spite's 25 years 
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ago, » 1.12 ± 0.38 x 1CT 10 as for example, 

7 Li/H = (1.23t{J;{g) x 1(T 10 (2) 

or 

7 Li/H = (1.1 - 1.5) x 1(T 10 

[22T |. Some uncertainty remains in the adopted stellar 
atmospheric temperatures, which may lead to somewhat 
higher estimates, as in the case of the globular cluster 
NGC 6397 

7 Li/H = 2.19 ±0.28 x 10~ 10 

[13, HI| or the somewhat controversial result 7 Li/H= 
(2.34±0.32) x 10" 10 ((H) for a sample of halo field stars. 
This is to be compared to the most recent standard BBN 
prediction of 

7 Li/H = (5.24t°; 74 ) x 10- 10 

[27| | taking into account experimental re-evaluations of 
the for the 7 Li abundance most important 3 He(a, 7) 7 Be 
reaction rate [H, [H, H3, [H|. Renewed measurements 
of this rate has not only increased the central value for 
7 Li/H by around 25% , but also reduced its error bar con- 
siderably, thereby leaving the discrepancy between pre- 
diction and observation more pronounced (around 4.2- 
5.3(7 [I?]] )■ It is therefore clear that an additional piece 
to the puzzle is missing. 

Whereas it seems essentially ruled out by now that 
either revised reaction rate data for 3 He(a,7) 7 Be or 
7 Be(Z?,p)2 4 He [32j or any other reaction, or a serious un- 
derestimate of the stellar atmospheric temperatures (of 
order 700K), or even a combination of the two effects 
may resolve the discrepancy, stellar depletion of 7 Li stays 
a viable, and non-exotic, possibility of resolving the pri- 
mordial 7 Li problem. Indeed, atmospheric 7 Li may be de- 
stroyed by nuclear burning when transported towards the 
interiour of the star. This is also observed in stars with 
large convective zones, but not in those relatively hotter 
(T ~ 6000K) radiation dominated main-sequence turn- 
off stars, which are/were believed, after two decades of 
research [33[ , to preserve at least a large fraction of their 
initial atmospheric 7 Li. Arguments presented in favor of 
non-depletion are (a) the essentially uniform observed 7 Li 
abundance over a wide range of different stellar temper- 
atures and metallicities, a pattern which is very difficult 
to achieve when significant depletion was at work, (b) 
the absence of star-to-star scatter in 7 Li abundances ex- 
pected to arise when, for example, depletion is induced 
due to rotational mixing of stars with different rotational 
velocities, and (c) the claimed presence of relatively large 
amounts of the much more fragile 6 Li isotope in some of 
those stars (cf. [34] for a review). Thus, taking Occam's 
razor, many observers believe in the 7 Li abundance on 
the Spite plateau to be the primordial one. Neverthe- 
less, the situation is clearly more complicated, since when 
gravitational settling ( "atomic diffusion" ) of heavier ele- 
ments is included, 7 Li is predicted to be depleted, albeit 



with a pattern which is not observed. In order to repro- 
duce the observations, another turbulent mixing process 
of unknown nature has to be "fine-tuned" to account for 
the observed pattern [3^, H|| . However, recent observa- 
tions of the observed patterns in not only 7 Li, but also 
Ca and Ti, in the globular cluster NGC 6397, lend in- 
deed some support to this idea, though being currently 
not statistically overly significant (2-3<r). The combina- 
tion of atomic diffusion and turbulent mixing has thus 
been claimed to be able to account for a 7 Li depletion 
of factor 1.9 |37| . a factor which could go a long way to 
solve the discrepancy. 

In contrast to 7 Li, standard BBN (SBBN) yields only 
negligible production of 6 Li/H- 10" 15 - 10" 14 [38]. This 
is mainly due to the only 6 Li producing reaction in SBBN 
D(a, 7) 6 Li being a quadrupole transition. 6 Li as observed 
in higher metallicity disk stars and in the Sun is believed 
to be due to cosmic ray nucleosynthesis, either result- 
ing from spallation reactions p,a+CNO^LiBeB or cos- 
mic ray fusion reactions a + a — >Li. Since the time- 
integrated action of standard supernovae producing cos- 
mic rays is measured by metallicity, it had been a sur- 
prise that three groups independently confirmed the ex- 
istence of 6 Li in the atmosphere of the [Z] « —2.2 low- 
metallicity star HD84937 ( 6 Li/ 7 Li« 0.052 ± 0.019) [H, 
a value too high to be comfortably explained by stan- 
dard cosmic ray nucleosynthesis. The surprise was even 
bigger when the long-term_pfoneering observational pro- 
gram of Asplund et al. [22| indicated the existence of 
relatively uniform 6 Li/ 7 Li ratios ~ 0.05 in about ten 
low-metallicity stars, reminiscent of a primordial plateau 
(metallicity-independcnt). Nevertheless, the situation is 
currently unclear. The presence of 6 Li in these hot stars 
has to be inferred by a minute asymmetry in the atmo- 
spheric absorption profile due to the blend of both, the 
6 Li and 7 Li isotopes. Any individual claimed 6 Li detec- 
tion is therefore only at the 2 — 4er statistical significance 
level. Based on observations (of a star which originally, 
however, was not claimed to have 6 Li) and complete hy- 
drodynamic 3D non-equilibrium simulations of stellar at- 
mospheres, Cayrel et al. [40( have recently asserted that 
line asymmetries due to convective motions in the at- 
mospheres could be easily misinterpreted as atmospheric 
6 Li, and that the claimed abundances are therefore spu- 
rious. In contrast, Asplund et al. [22j would even infer 
higher 6 Li/ 7 Li ratios, if they were to utilise their own 3D 
hydrodynamic non-local thermal equilibrium simulations 
of line profiles. Further analysis is clearly required. 

Concerning cosmic ray production of a putative 
6 Li/Hw 5 x 10~ 12 , a cosmic ray energy of lOOeV per 
interstellar nucleon is required, whereas standard super- 
novae generated cosmic rays may provide only 5 eV per 
nucleon at such low metallicities [Z]« —2.75 (of the star 
LP 815-43) [ilj. Such a 6 Li abundance requires therefore 
a very non-standard early cosmic ray burst, preferentially 
acting at higher redshift [42| , possibly connected to radio- 
loud quasars and the excess entropy in clusters of galaxies 
[43| , or to a significant fraction of baryons entering very 
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massive black holes |4l|. The requirements are extreme 
and no compelling candidate has been identified. Alter- 
natively, it may be that the 6 Li is produced in situ, due 
to fusion reactions in solar flares [441 ]. 

The anomaly in the Li abundance and/or the pur- 
ported one in the 6 Li abundance, could also be due 
to physics operating immediately during or shortly af- 
ter the BBN epoch, possibly connected to the dark 
matter. It has been long known that 6 Li is easily 
produced in abundance in non-thermal 3 H(a,n) 6 Li and 
3 He(a,p) 6 Li reactions during and after BBN without 
much disturbing the other light elements. Here ener- 
getic 3 H and 3 He may be produced by 4 He spallation- 
or photodisintegration- processes induced by the decay 
of relic particles 0, |4(| |47| . l48j or residual annihilation 
of dark matter particles J49|. A m x w 10 GeV neutralino 
of the WMAP density, and with substantial hadronic s- 
wave annihilation, as sometimes invoked to explain the 
DAMA/Libra anomaly [50j, would, for example, synthe- 
size a 6 Li abundance well in excess of that in HD84937. 
Similarly, the 7 Li abundance may be also affected by par- 
ticle decay (but less by annihilation). Early attempts 
to explain the 7 Li anomaly by 7 Li photodisintegration 
induced by the electromagnetic decay of a NLSP to a 
gravitino at r ~ 10 6 sec @, were subsequently shown to 
be in disaccord with either, a reasonable lower limit on 
the D/H ratio or an upper limit on the 3 He/D ratio [Hlj . 
However, it was shown that the injection of extra neu- 
trons between the mid and end ~ 100 — 3000 sec of the 
BBN epoch can affect a substantial reduction of the final 
7 Li abundance [13] . Moreover, if these neutrons are in- 
jected energetically, as is the case during the decay of a 
weak-scale mass particle, their 4 He spallation may at the 
same time be the source of an appreciable 6 Li abundance. 

This led to the realization that mf ~ 1 TeV supersym- 
metric staus decaying into w 100 GeV gravitinos dur- 
ing BBN may resolve two lithium problems at once [47l |. 
Subsequent more detailed calculations of stau freeze-out 
abundances, life times, and hadronic branching ratios 
within the constrained minimal supersymmetric standard 
model (CMSSM) |52| and 15311 confirmed these findings 
(cf. also to Ref. [54| and [55j| ~) . 

The present paper presents a re-/extended analysis of 
supersymmetric (SUSY) scenarios with gravitino LSP's 
which may solve one, or both, lithium anomalies. The 
reasons for such a reanalysis are multifold. 

1. Reference (52j as well as [53| analyzed scenarios 
only for large, electroweak scale gravitino masses. The 
present paper studies the case of light ~ 10 MeV-1 GeV 
gravitinos as well. Light gravitinos are typically emerging 
in gauge-mediated SUSY breaking scenarios (GMSB) . A 
consistent analysis of GMSB scenarios is performed. Fi- 
nally, for heavier electroweak scale gravitino masses, as 
predicted in the CMSSM, an improved scan of gravitino 
masses is presented. 

2. The adopted hadronic branching ratios in [52j] 
(based on Ref. [j]) are somewhat rough approximations 



only. For example, as shown in [12|, the hadronic branch- 
ing ratio for the decay of the lighter stau, fi — > Grqq, 
is larger when production of qq pairs by intermediate 
virtual photons, and not only Z-bosons, is taken into ac- 
count as well. 

3. The qq's have been noted to be injected at much 
lower energy [l2j] than taken in the approximative treat- 
ment of [52| . In particular, it has been argued that for ac- 
curate BBN calulations the effective energy chad should 
be taken as the mean in the qq invariant mass m qq . 

4. It will be shown in Section II, that taking the mean 
invariant mass m qq as proposed in 12j represents an ap- 
proximation as well. Indeed, a fully proper calulation re- 
quires knowledge of the energy spectrum of the injected 
nucleons. 

5. It has been recently pointed out that charged elec- 
troweak mass scale particles, such as the f , when present 
during the BBN epoch, may lead to interesting catalytic 
effects [56l. l57l [58|. In particular, catalytic enhancement 
of 6 Li production [H, [5!| has an important effect for 
stau life times t x > 3000 sec, whereas catalytic enhance- 



ment of 7 Li destruction [6Cj is less important for stau 
NLSPs. The impact of both effects on the t x w 10 3 sec 
lithium-solving parameter space has already been inves- 
tigated in Ref. (6l| (see also [HI, [63[ for other papers on 
catalysis). The present paper takes full account of the 
proposed catalytic processes. 

6. Similarly, Pospelov 64, 6f| has made the interesting 
suggestion that 9 Be may be produced due to catalytic 
effects. These effects are included in the present study 
and their significance is adressed. 

7. Finally, the re-evaluation of the main 7 Li producing 
rate [H Hi IM HI SH , as well as a shift in fl b h 2 [3, 
as indicated in the introduction, has increased the pre- 
dicted standard BBN 7 Li abundance. These changes are 
included in the present paper as well as recent changes 
in some catalytic rections rates (66|. 



II. DEPENDENCE OF RESULTS ON NUCLEON 
SPECTRUM 



In [12j it was argued that the quantity 



(m qq ) 



r to t(f) 



dm qq m qq 



dT(f -> rGqq) 
dm„« 



(3) 



should be taken as the "effective hadronic energy" 
£"Had(T — > i~Gqq) for the computation of BBN yields 
with decaying particles. It was noted that this is con- 
siderably lower for hadronic decays of ~TeV stau's than 
the adopted (rof — m^)/3 in 52]. Though this is a good 
comment, the situation is much more complicated. In 
contrast to electromagnetic decays, there is no obvious 
sensible definition of an "effective hadronic energy" . For 
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FIG. 1: D/H (upper) and 7 Li/H (lower) abundances com- 
puted for a number of selected models with rrvf ~ 1 TeV and 
fig/i 2 « 0.1 (see text for further detail) when (a) the invari- 
ant mass approximation is used (green triangles) and (b) the 
analysis is performed using a realistic nucleon energy spec- 
trum (black dots). The left (right) plateau show models with 
approximate decay times r ~ 10 3 sec (r ~ 300 — 600 sec), 
respectively. They correspond to models in the blue (red) 
regions of Fig. [5] respectively. 



example, ten qq fluxtubes of Eqm — 100 GeV produce 
about 13.5 neutrons and protons (of varying energy), 
whereas one fiuxtube of Ecm = 1 TeV produces only 3.2 
nucleons. Similarly, the injection often 10 GeV neutrons 
at temperature T — 30 keV lead to the production of 4.3 
thermalized D, 14.4 neutrons, and 4.0 nucleon number 
A = 3-nuclei. This needs to be compared to 0.5 D, 1.45 
neutrons and 0.7 A = 3-nuclei for one injected 100 GeV 
neutron. In both examples one started with the same 
initial "hadronic energy" , with the result being vastly 
different. The number ANf asc (neutrons, D, 3 H, 6 Li, ...) 
of thermalized nucleons/nuclei i produced per X particle 
decay at temperature T may be written in the following 



form 



AN? 



= B h 



E 



*E S P } 



had 



dm 



dNi 



-(E 3 ,T) (4) 



where Pf ad {Ej) = (dNj/dEj) is the distribution func- 
tion of nucleons j (protons, neutrons) with initial energy 
Ej injected into the plasma by the X decay, normalized 
such that 



dE .phad {Ej) 



(5) 



gives the mean number of injected nucleons j 
per hadronic decay. In the above the quantity 
dN t casc /dNj (Ej , T) gives the produced (and thermalized) 
nucleons/nuclei i per injected nucleon j injected at ini- 
tial energy Ej, and is the hadronic branching ra- 
tio. For the computation of AN£ asc a detailed anal- 
ysis of the thermalisation of an injected nucleon j, as 
well as the thermalisation of all produced secondary (and 
higher) generations of non-thermal nucleons and nuclei 
due to energetic nucleon inelastic- and spallation- scatter- 
ing pro cesses on thermal nucleons and nuclei is required 

( C f. 13). 

The invariant mass m q q gives the energy Eqm of 
the qq fiuxtube in the qq center-of-mass reference sys- 
tem. This quantity does not, however, give the en- 
ergy of the fiuxtube in the cosmic rest frame (stau rest 
frame). The latter is much larger. Nevertheless, approx- 
imating Pj lad (Ej) by the nucleon energy distribution in 
the center-of-mass frame of qq, resulting after fragmen- 
tation of a qq fiuxtube of energy (m qq -), is still some- 
what useful, since in certain temperature ranges, such as 
50keV< T< 80keV, the BBN yields mostly only depend 
on the number of nucleons injected, as long as they are 
somewhat energetic Eki n ^ 1 GeV. Whereas the energy of 
the created nucleons is not an invariant of the reference 
system, the number is. However, even this corresponds 
to an approximation. 

For a limited number of scenarios we have computed 
a realistic nucleon spectrum. This spectrum was then 
used as input to perform accurate BBN calculations. The 
scenarios studied are ~ 1 TeV staus decaying between 
t w 300 — 1200 sec into gravitinos, with the decay pro- 
ducing flh 2 k 0.1 in gravitinos. Fig. [1] shows the yields 
of D/H and 7 Li/H as obtained either from the realistic 
nucleon spectrum or from the approximation of a nucleon 
distribution resulting from a qq pair of energy (m q q) in 
its center-of-mass frame. It is seen that the realistic spec- 
trum gives larger 7 Li and lower Deuterium abundances 
by ^ 10 — 15% as compared to the {m qq -} approximation. 
It is cautioned however, that the relative error is only 
small for comparatively early decays. In what follows, 
when not stated explicitly otherwise, all figures shown 
employ this approximation. 

All BBN calculations assume Vt^h? = 0.02273 and em- 

tthe recent reevaluation of the 3 He(a,7) 7 Be rate by 
For details concerning the BBN calculations the 
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FIG. 2: Parameter space in the Qxh 2 B h versus r plane where 
decaying relic neutral particles have impact on the 7 Li and 
6 Li abundancies following Ref. [47j. Here Qx is the frac- 
tional contribution to the critical density the decaying particle 
would have today if it wouldn't have decayed, h is the Hub- 
ble constant in units of 100kms -1 Mpc -1 , B^ is the hadronic 
branching ratio, and r is the particle life time. The color cod- 
ing of the areas is as follows: (red) 7 Li/H< 2.5xl0~ 10 ; (green) 
0.015 < 6 Li/ 7 Li< 0.3; and (blue) both 7 Li/H< 2.5 x 10~ 10 and 
0.015 < 6 Li/ 7 Li< 0.3. Other (relevant) constraints on light el- 
ements are taken to be D/H< 5.3 x 10 -5 and the helium mass 
fraction Y p < 0.258, as indicated by a line. The black solid line 
labeled D/H> 4 x 10~ 5 shows how lithium-friendly parame- 
ter space is reduced when a less conservative limit on D/H is 
applied. The figure also shows by the grey and blue points, 
following essentially lines, the prediction for stau-NLSPs in 
the CMSSM assuming a gravitino mass of mg = 50 GeV and 
the prediction for neutralino NLSPs in a GMSB scenario as- 
suming a gravitino mass of mg = 100 MeV, respectively. To 
produce this figure, the hadronic decay of a 1 TeV particle 
into a quark-antiquark pair has been assumed. Results for 
other initial states and particle masses Mx may vary by a 
factor of a few such that the figure should be interpreted as 
indicative only. 

reader is referred to [67j (and Irilll ). Catalytic rates are 
taken from the recent papers |59j and [66j |. Details rele- 
vant to the CMSSM and GMSB supersymmetric models, 
NLSP freeze-out abundances, as well as the calculation 
of hadronic branching ratios will be presented elsewhere. 
This analysis was performed using several public numer- 
ical codes: SuSpect [68[ (a supersymmetric particle spec- 
trum calculator) , micrOMEGAs |69l| (a dark matter relic 
density calculator), CalcHEP [70| (an automatic matrix 
element generator), and PYTHIA [7l[ (a Monte-Carlo 
high-energy-physics event generator). 

III. A GUIDE TO THE MODEL BUILDER 

In this section we discuss, in a generic way, the relic de- 
caying particle parameter space that is interesting for the 
7 Li, 6 Li (and 9 Be) abundances. The discussion will then 
prove useful in the following sections for the identifica- 




FIG. 3: Parameter space in the Qxh 2 -Tx plane where charged 
relic decaying particles have impact on the cosmic 7 Li, 6 Li, 
and (possibly) 9 Be abundances without violating other light- 
element constraints. Here M x — 1 TeV, By t = 10 -4 and 
hadronic primaries of 1 TeV qq fluxtubes have been assumed. 
The color code is as described in caption Fig. [2] The two 
dashed lines delimit the band 2 x 10" 14 < 9 Be/H< 2 x 10" 13 , 
whereas the solid line indicates D/H= 4 x 10 . 



tion of scenarios in our particular model, supersymmetry 
with gravitino LSPs. Fig[2] shows the parameter space 
in relic decaying particle abundance, hadronic branching 
ratio Bh, and X-particle decay time which has impact 
on the cosmic lithium abundances while respecting all 
other abundance constraints. Note that results are only 
dependent on the product of fl x h 2 Bh, where fix is the 
fractional contribution of the A"-particle to the present 
critical density, if it wouldn't have decayed. (The fl x 
may be converted to the X-particle number to entropy 
density n x /s via n x M x /s = 3.6639 x lO^GeVl^/i 2 , 
where M x is the X-particle mass.) Evidently the results 
depend as well on the hadronic branching ratio of the X- 
particle decays. To produce the plots in Figs (Hand [3] we 
made the assumption of hadronic decays described by a 
M x = 1 TeV X-particle decaying into a qq pair occuring 
with probability B^. This does not correspond exactly to 
the situation of NLSPs decays into gravitinos G, i.e. neu- 
tralino x Gqq and stau f — > Grqq, such that, for our 
analysis, the figures should be used only as guides. How- 
ever, making different assumptions for particle mass and 
initial (hadronic) post-decay state make in many cases 
only changes of factors of order 2 — 3. 

It is seen that for early decay times within a " banana" - 
shaped region the 7 Li abundance may be significantly 
reduced ( 7 Li/H < 2.5 x 10~ 10 red area). Here the up- 
per envelope of the area is defined by the constraint 
D/H < 5.3 x 10~ 5 and the lower envelope by 7 Li/H 
>2.5 x lCT 10 . A value of D/H = 5.3 x 1(T 5 may be 
already uncomfortably large. Therefore, the accordingly 
labeled solid line in Fig[2] indicates the less conservative 
D/H< 4.0 x 10~ 5 limit. The green area indicates the 
region where a 6 Li/ 7 Li ratio of 0.015 < 6 Li/ 7 Li < 0.3 has 
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been synthesized. Here the upper end of the range re- 
quires already some post BBN stellar 6 Li depletion (rel- 
ative to 7 Li). The figure illustrates that for tx ^ 10 3 sec 
there exists plenty of parameter space which may produce 
an observationally important 6 Li abundance by hadronic 
particle decays. Moreover, as advocated in Rcf. (42|, a 
region around tx ~ 1000 sec and £lxh 2 Bh « 2 x 10~ 4 
may resolve both lithium anomalies at once. 

In Fig. [2] it has been implicitly assumed that the decay- 
ing particle is neutral. In case it is charged, catalytic ef- 
fects may have a strong impact on the 6 Li abundance [5t| . 
particularly for small Bh- Moreover, catalytic effects may 
also lead to the production of 9 Be [6J], thou gh it is cur- 
rently not clear if this indeed happens (cf. [66(). Note 
that 9 Be production is impossible with hadronic effects 
only. Fig[3]shows the 7 Li, 6 Li, and (possibly) 9 Be friendly 
areas for a charged massive particle decaying during the 
BBN era. Here a hadronic branching ratio Bh = 10~ 4 
and particle mass Mx = 1 TeV have been assumed. In 
order to compute 9 Be/H ratios we follow the assump- 
tions in Ref . [64j , keeping in mind the possibility of very 
significant modifications [66|]. It is noted that in FigfJ] 
an ordinate Vlxh 2 Bh has been chosen, whereas in FigO 
the ordinate flxh 2 is shown for one particular choice of 
Bh [z2|. It is intruiging to remark that when the de- 
caying particle is charged, not only may the 7 Li and 6 Li 
anomalies be solved, but for the same parameters an ob- 
servationally important 2 x 10~ 14 < 9 Be< 2 x 10~ 13 abun- 
dance [73| may be synthesized as well. This is seen by the 
dashed lines in Fig. [3] passing through the doubly pre- 
ferred region at flxh 2 Bh ~ 3 x 10 and tx ~ 1500 sec. 

Though the results of this section are essentially 
generic, and may be used as guide lines for construct- 
ing other particular "lithium-friendly" decaying particle 
scenarios, FigfJ] also shows predictions within two partic- 
ular setups in supersymmetric extensions of the standard 
model. The sequence of grey points gives the predic- 
tion of stau-NLSPs in the CMSSM under the assump- 
tion of a gravitino mass of rriQ — 50 GeV. Such scenarios 
come very close to the doubly-preferred blue region at 
tx ~ 1000 sec, which illustrates that (within the approx- 
imations) the 6 Li and Li observed abundances could be 
modified in an observationally favored way by stau de- 
cay, albeit for a stau mass m? ~ 1 TeV unaccessible to the 
LHC. This scenario will be analyzed once more in Section 
IV. It is noted here that though it seems somewhat incon- 
sistent to show f NLSP decay in Fig. [2j due to catalytic 
effects which are not taken into account in that figure, 
conclusions concerning the 7 Li and 7 Li and 6 Li preferred 
regions are hardly modified when catalytic effects are in- 
cluded (aside from 9 Be). The sequence of blue points 
gives the prediction of neutralino NLSPs in a GMSB 
model assuming itLq = 100 MeV. Such NLSPs, charac- 
terized by fairly large freeze-out densities and hadronic 
branching ratios, though they would produce too much 
6 Li when decaying late, could solve the 7 Li problem for 
decay times tx ~ 100 — 400 sec. This possibility is stud- 
ied in detail in Section V. 



CMSSM, tanp=10, A Q =0, m 3/2 E[10,500]GeV 




FIG. 4: Parameter space in the CMSSM mo — mi/2 plane 
where NLSP decay during BBN results in observationally fa- 
vored modifications of the 7 Li and/or 6 Li abundances with 
respect to standard BBN. Here tan/3 = 10 has been assumed. 
The color coding indicate: (red) 7 Li/H< 2.5 x 10~ 10 , (green) 
0.015 < 6 Li/ 7 Li< 0.15, and (blue) 7 Li/H< 2.5 x 10~ 10 and 
0.015 < 6 Li/ 7 Li< 0.66. Other constraints are as in Fig. [2] 
White areas are ruled out by BBN. The area labeled "LEP" is 
excluded by LEP accelerator constraints, whereas the area la- 
beled "EWSB" does not lead to electroweak symmetry break- 
ing. The line labeled " X NLSP" and "f NLSP" delineates the 
parameter space where a neutralino \ is the NLSP from that 
where the stau f is the NLSP. The three lines on the left 
hand side give contours for the lightest neutral Higgs mass, 
as labeled. 



IV. LITHIUM AND HEAVY LSP GRAVITINOS 

Weak scale mass "heavy" gravitinos are generally ex- 
pected in the constrained minimal supersymmetric model 
(CMSSM) , which will be utilised for the present analysis 
with the gravitino mass, > 10 GeV taken as a free 
parameter. Other parameters in the CMSSM are: mo, 
m±/2 - at the GUT-scale unified scalar, fermionic (gaug- 
ino) soft supersymmetry breaking masses, respectively, 
Aq - trilinear couplings in the scalar sector (assumed to 
be zero throughout), and tan/3 - the ratio of the vacuum 
expectation values of the two Higgs required in super- 
symmetry. It is known that the CMSSM with weak- 
scale LSP gravitinos may lead to cosmologically inter- 
esting changes in the primordial 6 Li and 7 Li abundances 
when the superpartner to the tau lepton, the stau, is the 
NLSP [13, [H3, ■ This may be seen given their abun- 
dances, which can be approximated for moderate tan (3 



7 




10 2 10 3 10 4 10 5 , 

T(s) 

FIG. 5: Gravitino relic density Q^h 2 versus NLSP (stau) life 
time r corresponding to the CMSSM models studied in Fig. [4] 
Color coding is as in Fig. [3] except that "pink" shows regions 
where 7 Li/H< 2.5 x 10~ 10 and 0.15 < 6 Li/ 7 Li< 0.66 account- 
ing for the possibility of 6 Li depletion. The cross-hatched 
region shows parameter space where 2 x 10 -14 < 9 Be/H< 
2 x 10 -13 , though the estimate is uncertain. Other lines show 
(as labeled) light element abundance constraints, the scan 
range in m.1/2 and mg, as well as the by WMAP inferred 
range for the dark matter density. 

by 

n*fc a « (2.2-4.4) xlO"^^)*, (6) 

their typical hadronic branching ratios -B/j « 1CP 4 — 10~ 2 , 
and life times in the range Tf ~ 10 2 — 10 5 sec as con- 
trolled by Eq.@ with k = 1, n — 4. Note that the range 
in the approximation Eq.fJfJ]) takes into account possible 
co-annihilation with the other sleptons (but does not ac- 
count for co-annihilation with a neutralino or other spe- 
cific effects at larger tan/3, such as mixings, higgs reso- 
nances, enhanced couplings, etc.). When compared to 
Fig. [2] and taking tan/3 = 10, one finds for instance 
that staus in the mass range m? «1- 1.5 TeV, having 
Bh ^ 10~ 3 and decaying into gravitinos with m,Q « 50 — 
200 GeV pass right through the preferred (blue) region. 
A somewhat heavier stau ~ 2 TeV, with a Bh ~ 2 x 10 -3 
and a Slfh 2 Bh « 1 x 10~ 3 , would be only 7 Li friendly 
as can be seen from Fig. [3J for w 100 GeV. The 
lithium friendly parameter space for the CMSSM with 
tan/3 = 10 is shown in Figs. [4] and [5] The calculations 
leading to Figs. Q] and \5\ include the improvements (1)- 
(3), (5), and (6) of the list given in the Introduction, but 
not improvement (4) . Treating each case with a realistic 



nucleon spectrum is numerically too expansive, and will 
be only done for a few models below. Note that Figs. 2] 
and [5] really show results of higher-than-two dimensional 
parameter space. Even at fixed tan f3 and Aq the results 
in the too — TO1/2 plane shown in Fig. [4] are for a vari- 
ety of gravitino masses. In particular, blue areas in the 
too — rnx/2 plane are covering up green and red areas 
which result for different choices of tuq. The figure, and 
other figures which follow, thus show where 6 Li (green), 
7 Li (red), and 6 Li and 7 Li (blue) friendly regions are ex- 
pected, in case m,Q has been approximately chosen. In 
contrast, in white regions for no ttiq within the adopted 
range may light element abundance constraints be met. 

Figs . |H and [5] may be directly compared to Figs. 1 and 
2 of [53] • It is seen that the doubly preferred (blue) region 
is narrower in [52| compared to the present study. This 
is simply due to a finer sampling of the gravitino mass 
in the present paper. Furthermore, models with larger 
stau life times, Tf > 5 x 10 3 sec, are due to catalytic 6 Li 
and (possibly) 9 Be overproduction now ruled out Ullil. 
Finally, it is noted that in the CI3/2I1 - tnlsp plane the 
doubly preferred (blue) region has hardly moved. This 
is somewhat surprising since for ~ 2 larger Bh (point 
2 in Introduction), as is the case, one would expect for 
the region to move a factor ~ 2 lower in Q,f h 2 (and thus 
^3/2 h 2 ) as to not overproduce D. Nevertheless, this effect 
is counter-balanced by a lower effective (m q q) (point 3 
in Introduction), implying less "distortion" of the light 
elements than initially envisioned. 

Scenarios where stau NLSPs decay at around t w 
1000 sec into gravitinos, thereby solving both lithium 
problems at once, have the added, and totally acciden- 
tal, benefit of coming tantalizingly close to producing all 
the dark matter in form of warm gravitinos during the 
decays. This is nicely illustrated in Fig. [5] where the blue 
area just overlaps from below with the WMAP strip, and 
models at higher S7 3 / 2 /i 2 are ruled out by Deuterium over- 
production, D/H> 5.3 x 10~ 5 . It is therefore interesting 
to see if this conclusions survives when a calculation of 
improved accuracy is performed. The results in Fig. O 
rely on approximating the " hadronic energy" by (m q q) of 
Eq.([3]). When properly calculating the hadronic energy 
release employing a realistic energy distribution of the 
nucleons leads to a milder D overproduction constraint 
than expected from the {m qq ) approximation as shown in 
section II, thus pushing the blue area more upwards into 
the WMAP strip. This is seen in Fig. [1] where results for 
CMSSM (tan/3 = 10) scenarios with n^h 2 close to 0.1 
are shown, which have been computed without the m q q 
approximation. Such scenarios are thus fully consistent 
with producing all the dark matter non-thermally. 

It is of interest if catalytic effects due to the electri- 
cally charged staus may also lead to a cosmologically 
important 9 Be abundance [H, [65[ in the doubly pre- 
ferred parameter space. Indeed, this is the case, as may 
be seen in Fig. [5] where the 9 Be/H -ratio in the range 
2 x 10~ 14 < 9 Be/H< 2 x 10~ 13 is indicated by the cross- 
hatched region. This is indeed very interesting since 
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the observed 9 Be/H in the lowest metallicity stars is 
9 Be/H« 3 x lCT 14 - 1(T 13 [73], thus close to the pre- 
dicted one. 

In Table I parameters and BBN yields for some partic- 
ular SUSY points in the CMSSM (and GMSB, cf. Sec- 
tion V) for tan/3 = 10 are shown. One notes two classes 
of models, those which synthesize much 6 Li (and 9 Be), 
for larger tx, and those which significantly reduce 7 Li 
for smaller tx, with the latter models also accounting 
for the dark matter produced during the decay. In ei- 
ther model D/H may be significantly less than 4 x 10~ 5 . 
Models which accomplish both, a factor > 2 destruction 
of 7 Li as well as cosmologically interesting production of 
6 Li (and possibly 9 Be), however, typically tend to have 
D/H> 4 x 1CT 5 . Finally, it is noted that even in the pres- 
ence of hadronic decays, i.e. injection of extra neutrons, 
the 9 Be/ 6 Li ratio falls typically in the range 1 — 4 x 1CP 3 , 
as proposed in Ref. [64| in the absence of extra neutrons. 

Last but not least Fig. [5] shows results for the CMSSM 
and tan/5 = 50. As was already noted in (52j . for tan/3 
large, and a stau NLSP (for neutralino NLSPs cf. Section 
V), typical flf h 2 Bh(f) < 10 -4 are too small to resolve 
the 7 Li problem (cf. Fig. This is due to the smaller 
freeze-out abundance due to efficient annihilations of f f 
into W + W~(ZZ) and hh (light CP- even Higgses) for 
large tan/3, as compared to the ff — > rr channel which 
controls the abundance at low tan/3. Such effects can 
be due to Higgs poles or large left-right stau mixing for 
relatively light stau masses (as shown in [76|, [77| in non- 
CMSSM scenarios). In the CMSSM and for very heavy 
stau NLSPs as is the case in Fig. d W+W~ plus ZZ 
channels become dominant. Nevertheless, there remains 
parameter space which may synthesize observationally 
important primordial 6 Li and 9 Be abundances at once. 
This may occur for staus as light as rrif « 500GeV, po- 
tentially visible at the LHC. It is noteworthy that for 
very heavy rrif > 2TeV there are again regions where ei- 
ther 7 Li alone or both 6 Li and 7 Li are synthesized at the 
observationally inferred level. Furthermore, due to the 
smallness of ftf h 2 for relatively light staus, only a small 
fraction of the dark matter would be produced by stau 
decays, with the "missing" gravitinos possibly produced 
during reheating at a comfortably large reheat tempera- 
ture of Tjijj ~ 10 10 GeV. Very heavy stau decays can still 
account for all the (gravitino) dark matter as shown in 
Fig. nevertheless models at smaller tan (3 ~ 10 seem 
more economical in the context of dark matter genera- 
tion. 



CMSSM, tan|5=50, A =0, m s=[10,500]QeV 




10- 



10 J 



m 1/2 (GeV) 
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FIG. 6: Same as Figs|4]and[5] but for tan/3 = 50. The region 
labeled "tachyonic" is excluded due to some sfermion masses 
becoming negative. 



V. LITHIUM AND LIGHT GRAVITINOS 

In this section we identify SUSY parameter space 
which results in significant changes of the primordial 
6 Li and 7 Li abundances, in cases when the gravitino is 
rather light < 10 GeV. Such an analysis has so far not 
been performed. Light gravitinos are a typical prediction 
of gauge- mediated SUSY breaking scenarios (GMSB), 



where SUSY breaking in a hidden sector is transmitted 
to the visible sector via gauge interactions of some mes- 
senger fields [zll , [7!| . Assuming a typical grand unified 
group, such models have three continuous and one dis- 
crete parameter, namely A, the common mass scale of 
the soft Susy breaking masses, M mess , the Susy preserv- 
ing messenger mass, tan /3, the ratio of the two higgs 
doublet vevs and N mess the number of quark-like (as- 
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TABLE I: Light-element abundances yields and gravitino 
abundance Q, 3 / 2 h 2 for a number of selected models. Cf. to 
Table II for the particular particle model parameters. Model 
denoted by (a) have been computed with the invariant mass 
approximation (see text for details). 



Model 


D/H 


7 Li/H 


6 Li/ 7 Li 


9 Be/H 


«3/2^ 2 


A 


3.45 ■ 10" 


■ r > 


2.20 • 10" 


10 






0.11 


A" 


3.86 ■ 10" 


5 


1.66- 10" 


10 






0.11 


B 


4.75 • 10" 


5 


2.14 • 10" 


10 


0.058 


3.4 ■ 10~ 14 


0.10 


B a 


5.42 • 10" 


5 


1.83- 10" 


10 


0.081 


3.8 ■ 10~ 14 


0.10 


C a 


3.59 ■ 10" 


5 


3.32 • 10" 


10 


0.044 


7.5 • HT 14 


0.0086 


D a 


2.57- 10" 


5 


4.95 ■ 10" 


10 


0.044 


1.9 • 10 -13 


3- 10~ 4 


E a 


3.61 • 10" 


5 


2.09 • 10" 


10 






2 • 10 -4 



TABLE II: Particle physics model parameters corresponding 
to the abundance yields shown in Table I. All masses are in 
GeV and tan (3 = 10 



Model 


SUSY - NLSP 


Mnlsp 


m G 


TX (s) 


A 


CMSSM f 


2060 


160 


403 




(m ,m 1/2 ) = (1045,5020) 








B 


CMSSM f 


1638 


175 


1553 




(m ,m 1/2 ) = (678,4200) 








C 


CMSSM f 


763 


30 


1980 




(m ,m 1/2 ) = (395,1831) 








D 


GMSB N = 2f 


264 


2.72 


3320 




(A,M mcss ) = (10 5 ,5- 10 6 ) 








E 


GMSB N = 1 x 


133.8 


0.08 


103 




(A,M mcss ) = (10 5 ,5- 10 6 ) 









sumed to be equal the number of lepton-like) multiplcts 
(of some GUT group) messenger supermultiplcts. The 
soft Susy breaking gaugino (square of scalar) masses are 
then generated at the 1- (2-) loop level, respectively, lead- 
ing to the physical spectra and couplings of the MSSM 
particles. (The soft breaking trilinear couplings A are 
generated only at the two-loop level and will be set to 
zero). Even though the fine details of the MSSM spec- 
trum and couplings depend on these four parameters, in 
most parts of the parameter space varying M mcss , which 
affects the soft masses only logarithmically, will be effec- 
tively irrelevant to our study. We are thus left with one 
free mass parameter A which should furthermore be of 
order 100 TeV if the MSSM spectrum is to remain at the 
electroweak mass scale. Such models are therefore more 
constrained than, for example, the CMSSM which has 
(at least) two mass parameters mo and TO1/2 (given that 
trilinear couplings A's have been set to zero). Neverthe- 
less, before discussing results in the GMSB, we will still 
study some aspects of the CMSSM when the gravitino 



mass is (arbritrarily) low. Though not quite consistent, 
as one expects m^ ~ rn so f t when SUSY breaking in a 
hidden sector is communicated to the visible sector by 
gravitational interactions, the study of the CMSSM with 
light gravitinos may help to shed more light on existing 
" lithium- friendly" parameter space in models beyond the 
toy models CMSSM and GMSB. 

Consulting Fig. [2] as the key figure of where to ex- 
pect parameter space solving the 7 Li anomaly, one ob- 
serves that for all 3 x 10~ 4 < ^NhSph 2 Bh ^ 0.1, and with 
a somewhat tuned gravitino mass to match the " desired" 
NLSP life time, solutions should exist. It is noted here 
that the exact range required in ^nlsp^-B/i depends on 
the NLSP mass Mnlsp, moving down as Mnlsp moves 
down (Fig. [2] assumes Mnlsp = 1 TeV) and vice versa. 
As discussed above, stau NLSP abundances are not much 
larger than fif h 2 ~ 1 (cf. Eq. 0, whereas their hadronic 
branching ratios are usually small Bh < 10 -3 . Thus, 7 Li 
solving areas for the stau maybe only found for large 
TOf ~ 1 TeV since only there ttfh 2 ~ 1. Furthermore, 
to match the desired life time, large mf require large 
nig (cf. Eq. ((21)), thus pointing towards gravity medi- 
ated susy breaking scenarios. These solutions have been 
all discussed in Section IV. In contrast, the freeze-out 
Q. x h 2 for neutralinos may be large ~ 1 — 100 and there 
hadronic branching ratio due to x — > Gqq with an inter- 
mediate Z or photon is typically of the order of Bh ~ 0.1. 
In particular, one may find 7 Li-friendly models for com- 
paratively light neutralinos only for which fl x h 2 and Bh 
are not too large as to saturate the upper bound on 
^NLSpBh(NLSP) given by 4 He overproduction. Since 
t x scales with m x to the fifth power, the gravitino mass 
m,Q in such decays has to be rather light in order to match 
the decay time window 100 sec < r < 1000 sec (cf. Fig. [2]). 
Thus, one is automatically led to lighter gravitinos. This 
implies also, since Qnlsp may not be too large, and 
^3/2 = ^NLSp{ m Ql m NLSp) that only a small fraction 
of the dark matter would be created by such decays. 

Identifying areas of small neutralino abundance 
in the CMSSM in the mo-m 1 / 2 plane, in or- 
der to satisfy n NLS pB h (NLSP) < 0.1 (rather 
Q N LSpB h (NLSP)< 0.03 for m x w 100 GeV), is remi- 
nescent of identifying the regions where fl x ~ Qwmap 
to account for neutralino dark matter. The latter 
question has received much attention over the last 
years. It is therefore not surprising that one finds 
potential 7 Li-solving parameter space in three distinct 
and "known" regions (a) the bulk region close to the 
LEP bound at small m 1/2 « 300 - 400 GeV, (b) the 
co-annihilation region close to the line of neutralino-stau 
mass degeneracy, and (c) the focus point region at 
large mo and m 1 / 2 - These regions are shown in red 
in Fig. [7] for tan/3 = 10. In contrast with the heavy 
gravitino case, here Higgs mass bounds from LEP could 
potentially exclude large fractions of the bulk region, as 
illustrated in the figure with lines of constant lightest 
Higgs mass including the Standard Model Higgs present 
lower limit m^o ^ 114 GeV. Although the experimental 
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CMSSM, tanp=10, A Q =0, m 3/2 s[0.01,10]GeV 




FIG. 7: 7 Li (red) and 6 Li (green) solving regions for a 
CMSSM model at tan/3 = 10 and with light graviti- 
nos in the range 10 MeV ^ ^ 10 GeV. All constraints 
and lines have similiar meaning to those in Fig. [4] Note 
that the 7 Li solving regions exist essentially whereever 3 x 
HT 4 < fi N LSpk 2 Bh(NLSP) < 0.1. 




FIG. 8: Parameter space in the gravitino mass nig NLSP 
mass tonlsp plane for light gravitinos < 10 GeV which 
significantly reduce the 7 Li abundance, i.e. 7 Li/H< 2.5 x 
10 -10 (red), while respecting all other BBN constraints, as 
labeled. Here a GMSB model with tan/3 = 10, N mess = 1, 
and M m esB = 5 x 10 6 GeV has been assumed. Light-blue areas 
are allowed by BBN while white areas are ruled out. 



limits from LEP and Tevatron are milder for the MSSM 
higgses, one should keep in mind that a detailed study 
is needed, beyond the maximal mixing or no-mixing 
assumptions, in order to assess quantitatively the effect 
of experimental exclusions on our scenarios. There 
exists more parameter space to solve the 7 Li problem by 
neutralino decay than to produce neutralino dark matter 
of the right density (when uIq > m x ). Here region (a) 
requires gravitino masses 30MeV < m<=j < lOOMeV, region 
(b) 200MeV < m^ < 6.5GeV, whereas region (c) prefers 
200MeV< m 6 < 800MeV. Note also that for larger tan/3 
values ( > 40) a fourth 7 Li friendly funnel shaped region 
appears, corresponding to the well-known CP-odd higgs 
resonance effects. 

We next study GMSB models. The blue sequence 
of points in Fig[2] shows the prediction for VL x h 2 Bh for 
NLSP ncutralinos decaying into 100 MeV gravitinos in a 
particular GMSB model, employing N mess = 1 (lepton- 
and quark-like) messenger particles of mass M mess = 
5 x 10 6 GeV and tan/3 = 10. More accurate results for 
this model, now with varying mg, are shown in Fig. [8] 
There clearly exists 7 Li-solving parameter space for light 
neutralinos < 200GeV detectable at the LHC. We have 
also shown in light blue the parameter space region con- 
sistent with standard BBN (whereas the white region is 
ruled out). This is interesting by itself, as it illustrates 
the overall consistency of GMSB scenarios with primor- 
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FIG. 9: As Fig. [8] but with N mesa = 2, leading to neu- 
tralino NLSPs for m N LSP ~ 175 GeV and stau NLSPs for 
tinlsp ~ 175 GeV as indicated. The green area shows regions 



with 0.015 < 6 Li/ 7 Li < 0.66. 
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TABLE III: Potential for SUSY with gravitino LSPs to re- 
solve the 7 Li problem, account for 6 Li, produce 9 Be, account 
completely for the dark matter due to non-thermal decay pro- 
duction, and be detectable at the LHC. 



Gravitino 


NLSP 


7 Li 


6 Li 


9 Be 




LHC 


light 


stau 
neutralino 


X 

V 


V 
X 


V? 
X 


X 
X 


V 
V 


heavy 


stau 
neutralino 


V 

X 


V 

X 


V? 
X 


v 7 
X 


X 

V 



dial nucleosynthesis irrespective of whether the solution 
to the 7 Li problem is of particle physics or other as- 
trophysical origin. In particular a combination of the 
standard bounds on 4 He mass fraction and Deuterium 
abundance cuts all gravitino masses > 200 MeV for a neu- 
tralino NLSP in the "natural" range < 1 TeV, thus rein- 
forcing the common lore of light gravitinos within GMSB. 
Finally, we illustrate in Fig. [9] a GMSB scenario with two 
messenger (super)multiplets, N mess = 2. In this case, the 
NLSP can be a stau or a neutralino, depending on the 
values of A. Here a 6 Li-solving pattern requires a stau 
NLSP heavier than ~ 175 GeV with m & > 1 GeV, while 
solutions for 7 Li still obtain for a neutralino NSLP in the 
same part of the parameter space as in the N mess = 1 
case. The discontinuity in Fig. [5] is obviously due to the 
different nature of the NLSP, different hadronic branch- 
ing ratios as well as absence of catalysis effects for neu- 
tralino NLSP. We note also that an intermediate range 
of gravitino masses is disfavoured if one insists on solv- 
ing at least one of the two lithium problems on top of 
the consistency with standard BBN. However this gap in 
the mass range is obtained for a fixed messenger mass. 
We checked that it can be easily filled by varying the lat- 
ter over a few orders of magnitude. Other patterns can 
appear by varying the number of messengers, or by re- 
laxing the grand unification assumptions taking different 
numbers of quark-like and lepton-like messenger fields. 

VI. CONCLUSIONS 

The present paper has studied a variety of super- 
symmetric scenarios with gravitino LSPs. It has been 
mostly motivated by an apparent mismatch between the 
by standard BBN (SBBN) predicted 7 Li/H ratio and 
that observed in low-metallicity halo stars, as well as the 
(potentially controversial) claims of an observed 6 Li/ 7 Li 
plateau at low metallicity reminiscent of a 6 Li primor- 
dial abundance factor ~ 10 3 larger than predicted by 
SBBN. Though one or both of these anomalies may have 
an astrophysical/observational origin, it is conceivable 



that they point to physics beyond the standard model 
of particle physics and BBN. In particular, supersym- 
metric (NLSP) particle decays during BBN may explain 
either one or both of these anomalies. Though not the 
first such study, the present study improves in accuracy 
on several accounts, as well as studies, for the first time, 
the case of light gravitinos 10 MeV < m^ < lOGeV. De- 
tailed and improved analyses of relic NLSP abundances 
and hadronic branching ratios are performed in two dis- 
tinct SUSY models, the CMSSM with heavy gravitinos 
m 6 >10GeV (Section IV) and the GMSB with light 
gravitinos 10 GeV (Section V). Other significant 

improvements concern the study of BBN yields due to 
NLSP hadronic decays when no simplifying assumptions 
are made about the injected nucleon spectrum (Section 
II), as well as the full inclusion of catalytic effects dur- 
ing BBN due to the presence of weak-mass scale electri- 
cally charged particles (i.e. the stau). Finally, approxi- 
mate "model builders instructions" of where, on general 
grounds, to find relic particle decay parameter space rel- 
evant to the 7 Li and 6 Li abundances are given as well 
(Section III). 

The general results of this study are summarized in 
Table III. Susy with either heavy or light gravitino LSPs 
may impact the cosmic lithium abundances while satis- 
fying all other BBN constraints. Scenarios with heavy 
gravitinos require the charged stau to be the NLSP, and 
may solve the 7 Li problem with or without production 
of Li in abundance as claimed to be observed in low- 
metallicity stars, while producing all the dark matter 
as inferred by WMAP. These conclusions are unchanged 
from earlier ones, independent of catalysis effects and 
improvements in hadronic branching ratios and injected 
nucleon spectrum due to the decay. However, in those 
parts of the parameter space where significant 6 Li is syn- 
thesized, interestingly, it is conceivable that 9 Be on levels 
consistent with observations is produced due to catalytic 
effects as well. Unfortunately, scenarios of stau NLSP 
with a heavy gravitino LSP are most likely untestable at 
the LHC, since m? ^ 1 TeV is preferred. This is different 
for scenarios with light gravitinos. Here, either cosmo- 
logically important abundances of 6 Li (and possibly 9 Be) 
may be synthesized (for stau NLSP) or the 7 Li problem 
may be solved (for neutralino NLSP) for NLSPs light 
enough to be detectable at the LHC. However, in such 
scenarios only a small fraction of the dark matter is pro- 
duced, with the remaining gravitino dark matter possibly 
produced during a reheating period after inflation. 
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